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ABSTRACT 

The  effect  of  cornercube  retroreflector  orientation  on  satellite  laser  ranging  (SLR)  detected 
photon  counts  is  examined.  It  is  found  that  a  retroreflector  tilt  can  ameliorate  the  effects 
of  the  satellite  velocity  aberration  by  broadening  the  reflected  beam  through  diffraction. 
The  optimal  tilt  angle,  away  from  the  SLR  station,  can  be  up  to  approximately  20.5°  for 
the  present  retroreflector  design.  The  size  and  direction  of  the  tilt  angle  depend  on  the 
satellite  velocity  relative  to  the  SLR  station,  and  the  direction  of  the  satellite  from  the  SLR 
station.  The  required  attitude  control  tolerance  is  to  within  17°  of  the  optimal  attitude 
control  strategy  determined  in  the  present  work.  A  pre-launch  measurement  of  the  re¬ 
flectance  (diffraction)  pattern  of  each  retroreflector  is  recommended  to  take  account  of 
the  significant  likely  manufacturing  variation. 
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Satellite  Laser  Ranging  Photon-Budget  Calculations 
for  a  Single  Satellite  Cornercube  Retroreflector: 
Attitude  Control  Tolerance 


Executive  Summary 

An  early  design  study  [1]  in  2011  considered  the  suitability  of  a  corner-cube  retroreflector 
of  radius  7.5  mm  for  a  small-satellite  mission.  Some  estimates  of  the  laser  beam  transmit¬ 
ter  and  receiver  optics  were  assumed,  and  pointing  angles  from  the  Zenith  to  60°  were 
considered.  The  retroreflector  was  found  to  give  an  adequate  number  of  detected  pho¬ 
tons  for  a  reasonable  range  of  atmospheric  conditions. 

The  present  study  concerns  a  smaller  solid  corner-cube  retroreflector  (with  a  radius  of 
6.35  mm)  that  has  been  adopted  in  the  current  satellite  design.  In  particular,  the  study 
has  concentrated  on  determining  the  tolerances  on  the  retroreflector  orientation. 

The  laser-beam  transmitter  and  receiver  optics  parameter  values  have  also  been  refined, 
as  have  the  tolerances  on  pointing  accuracy  and  these  have  been  incorporated  in  the 
present  study.  These  include  a  smaller  beam  divergence  and  greater  pointing  accuracy, 
but  lower  detector  efficiency  and  optics  transmittance  and  greater  range  of  pointing  an¬ 
gles  (up  to  75°  from  the  zenith). 

Two  important  aspects  of  the  dependence  of  the  detected  photon  count  on  the  retroreflec¬ 
tor  orientation  have  been  included  in  the  present  study  that  were  omitted  from  the  earlier 
work.  The  first  is  the  shift  in  reflected  beam  angle  due  to  the  satellite  speed;  an  effect 
known  as  the  velocity  aberration.  In  the  present  case,  this  causes  a  shift  in  beam  direction 
of  up  to  11  arcsec.  The  size  of  the  aberration  depends  on  the  angle  between  the  satel¬ 
lite  velocity  and  the  direction  from  the  satellite  to  the  range  station,  and  so  varies  with 
pointing  angle.  The  second  is  the  orientation  dependence  of  the  retroreflector  diffraction 
pattern.  This  is  a  non-trivial  exercise  due  to  the  complications  involved  in  calculating 
the  apparent  aperture  size  of  the  retroreflector  aperture.  Most  of  the  available  literature 
concerns  retroreflectors  with  a  minimal  depth-to-radius  ratio  (of  1.41),  whereas  the  cur¬ 
rent  retroreflector  of  interest  has  a  larger  ratio  (1.85)  and  so  the  theoretical  basis  of  the 
modelling  has  been  re-derived. 

It  is  found  that  the  retroreflector  should  be  tilted  away  from  the  SLR  station  in  order  to 
broaden  the  reflected  beam  by  diffraction  from  a  narrower  aperture.  This  mitigates  the  ef¬ 
fects  of  the  velocity  aberration.  The  optimum  tilt  angle  varies  with  the  size  of  the  velocity 
aberration.  At  the  maximum  expected  velocity  aberration  of  10.92  arcsec,  a  retroreflector 
tilt  of  (20.5  ±  17)°  is  required  to  exceed  the  unity  photon-count  threshold,  with  a  tilt  of 
20.5°  being  optimal.  The  axis  of  the  rotation  is  perpendicular  to  both  the  satellite  (relative) 
velocity  and  the  direction  to  the  SLR  station. 
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The  same  tilt  angle  (and  tolerance)  can  be  used  at  lower  velocity  aberration  values  and 
still  exceed  the  detection  threshold.  Nevertheless,  a  reduced  tilt  angle  (according  to  the 
values  in  Table  2  and  plotted  in  Figure  9)  will  give  optimal  photon  counts. 

The  manufacturing  tolerances  of  the  cornercube  retroreflectors  identified  in  the  current 
satellite  design  are  found  to  allow  for  a  significant  variation  in  the  reflected  beam  width. 
The  possible  diffraction  patterns  encompass  a  range  of  different  attitude  control  strate¬ 
gies.  It  is  likely  that  a  different  strategy  will  be  required  for  each  cornercube,  depending 
on  its  particular  diffraction  pattern.  It  is  even  possible  that  some  retroreflectors  might 
not  be  suitable  for  use  in  the  current  application.  It  is  therefore  recommended  that  the 
pattern  of  each  retroreflector  be  measured  and  an  attitude  control  strategy  developed  to 
match  each  one. 
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1  Introduction 

The  Biarri  satellite  will  have  a  cornercube  retroreflector  mounted  in  its  surface.  This  is  to 
enable  accurate  ranging  experiments  to  be  performed  on  it  using  satellite  laser  ranging 
(SLR)  techniques.  The  SLR  will  be  conducted  from  Mount  Stromlo  in  eastern  Australia. 
An  early  design  study  [1],  hereafter  referred  to  as  S&P,  found  that  a  single  cornercube 
retroreflector  with  aperture  radius  of  7.5  mm  would  give  a  sufficient  reflected  signal  for 
SLR  using  the  transmission  and  receiving  SLR  equipment  at  Mount  Stromlo. 

The  aim  of  the  present  work  is  to  determine  the  attitude  control  necessary  to  ensure  an 
adequate  signal  from  the  current  retroreflector  planned  for  the  satellite.  The  signal  is  re¬ 
duced  if  the  retroreflector  aperture  is  tilted  away  from  the  SLR  station  (hereafter  referred 
to  simply  as  the  station)  due  to  the  smaller  projected  aperture  area.  There  are  other  effects 
that  also  need  to  be  considered.  The  diffraction  pattern  of  the  reflected  light  changes  with 
retroreflector  tilt,  and  an  effect  known  as  the  velocity  aberration  will  rotate  the  reflected 
distribution  of  light  due  to  the  satellite's  high  speed.  These  two  effects  are  considered 
in  some  detail  in  the  present  work.  It  is  found  that  there  is  an  optimal,  non-zero  tilt  of 
the  retroreflector  away  from  the  direction  to  the  station.  This  tilt  angle  depends  on  the 
satellite  location  (relative  to  the  station),  and  the  attitude  tolerances  are  quite  wide  (of  the 
order  of  10°)  if  the  detection  threshold  is  to  be  met  or  exceeded. 

There  are  two  notable  effects  that  have  not  been  included  in  the  present  study.  These  are 
the  effects  of  atmospheric  turbulence  on  the  propagation  of  the  laser  beam  on  both  out- 
and  in-bound  journeys,  and  the  sensitivity  of  the  reflectivity  of  the  cornercube  faces  to 
different  states  of  the  laser  polarisation.  The  former  will  cause  fluctuations  in  the  signal 
strength,  but  are  expected  to  preserve  the  average  values  calculated  in  the  present  work. 
The  latter  might  be  significant  in  some  circumstances  of  retroreflector  attitude  where  the 
reflection  angles  are  high  and  there  is  an  approximately  10%  variation  on  reflectance 
between  extreme  states  of  polarisation  [2],  These  situations  are  expected  to  be  transitory, 
and  so  not  of  significance  in  these  photon-budget  calculations  where  an  average  of  the 
two  reflectance  values  has  been  used. 


2  SLR  geometry  and  detection  threshold 

The  geometry  of  the  components  of  an  SLR  operation  is  illustrated  in  Figure  1.  At  the 
Mount  Stromlo  station  the  transmitter  and  receiver  are  co-located  (sharing  the  same  fore¬ 
optic)  so  that  the  receiver  angle  0r  =  ±a,  where  the  sign  of  the  velocity  aberration,  a, 
depends  on  the  direction  of  the  satellite  velocity  (as  discussed  in  §4.2.1).  A  laser  beam 
pulse  is  sent  from  the  transmitter  toward  the  retroreflector  and  the  time  taken  for  the  first 
return  photon  to  be  detected  by  the  receiver  is  used  to  determine  the  distance  travelled, 
and  hence  the  range  of  the  retroreflector.  The  pointing  angle  <i>  <  75°,  and  0,  is  the  retrore¬ 
flector  tilt  angle;  the  parameter  of  primary  interest  in  the  present  work.  The  remaining 
parameters  are  defined  and  discussed  at  greater  length  in  the  following  sections. 

The  details  of  how  the  pulse  travel  time  is  measured  precisely  is  not  of  concern  in  the 
present  work,  which  is  focused  mainly  on  modelling  the  strength  of  the  return  signal. 
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Figure  1:  Diagram  (not  to  scale)  showing  the  angles  defined  for  the  satellite  laser  ranging 
geometry.  The  transmitter  (Tx)  and  receiver  (Rx)  are  in  general  separate  (as  shown),  but  in 
the  present  application  will  be  co-located.  The  satellite — and  hence  the  retroreflector — is 
moving  with  velocity  v  relative  to  the  ground  station.  The  velocity  component  normal  to 
the  displacement  vector  r  between  the  transmitter  and  retroreflector  is  v' .  The  reflected 
light  distribution  is  rotated  through  angle  a  by  the  effect  known  as  the  velocity  aberration. 
If  the  transmitter  and  receiver  are  co-located  then  the  receiver  angle,  9r  =  —a. 


Nevertheless,  it  is  worth  noting  that  the  detection  of  the  return  signal,  and  distinguish¬ 
ing  it  from  background  noise,  is  a  statistical  process  that  requires  multiple  pulses.  The 
required  minimum  detectable  photon  count  per  pulse  [3]  is  approximately  np  >  0.02  (2 
photons  per  100  pulses)  at  night  and  np  >  0.5  (1  photon  per  2  pulses)  during  daylight. 
For  the  current  calculations  a  unity  threshold  (np  =  1  photon  per  pulse)  is  therefore  as¬ 
sumed,  both  for  simplicity  and  to  allow  a  margin  of  error  for  the  many  uncertain  variables 
involved,  such  as  the  atmosphere. 


3  Retroreflector  design 

The  retroreflector  planned  to  be  used  on  the  satellite  is  a  solid,  crown-glass  cornercube 
with  silvered  reflecting  surfaces.  It  is  an  off-the-shelf  item  from  a  large,  optical-component 
manufacturing  company.  The  retroreflector  design  parameter  values  are  listed  in  Fig¬ 
ure  2.  The  crown-glass  material  N-BK7  has  a  refractive  index  [2]  n  =  1.519  at  the  wave¬ 
length  (A  =  532  nm)  used  for  SLR  from  Mount  Stromlo.  The  surface  silvering  improves 
the  reflectance  over  uncoated  glass  surfaces.  It  is  assumed  in  the  present  study  that  the 
reflectance  of  the  silver  coating  will  be  equal  to  that  of  silver  metal  (0.95)  [2],  giving  a  nett 
three-surface  reflectance  for  the  cornercube  of  p  =  0.857.  The  circular  aperture  (radius 
rcc  =  6.35  mm)  of  the  retroreflector  is  not  recessed  in  its  holder  or  the  satellite.  This  is 
important  in  determining  its  reflective  area  at  oblique  ranging  angles,  and  its  consequent 
reflected  diffraction  pattern. 
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Edmund  Optics  #45-203 

Outer  diameter,  OD' 

31.75  mm 

Inner  diameter,  ID  (=  2rcc) 

12.70  mm 

Height,  H  (=  l) 

11.73  mm 

Length,  X+ 

18.08  mm 

Diameter  tolerance 

±0/  -  1  mm 

Surface  accuracy 

1/8  A 

Surface  quality 

60-40 

Housing  tolerance 

ODd  ±0/  —  0.5  mm 
H:  ±0.25  mm 

Beam-angle  tolerance 

3  arcsec 

Substrate 

N-BK7 

Coating 

Internal  silver 

Figure  2:  Design  diagram,  representative  photographs,  and  a  list  of  physical  properties  of 
the  retroreflector  solid  cornercube  under  consideration  (from  [4]).  The  quantities  marked 
"1"  do  not  directly  affect  optical  performance. 


3.1  Cornercube  retroreflector  dihedral  angle  tolerance 

The  reflected  diffraction  pattern  from  a  cornercube  retroreflector  depends  critically  on 
the  sizes  of  the  three  dihedral  angles  between  the  reflecting  faces.  If  these  angles  are 
exactly  right  angles,  then  the  narrow-beam  (that  is,  without  the  diffraction-broadening 
effects)  reflected  direction  will  be  exactly  the  opposite  of  the  incoming  beam,  as  required. 
Offsets  in  these  dihedral  angles  will  in  general  lead  to  a  deviation  from  the  ideal  return 
beam  direction.  As  is  often  the  case,  the  present  cornercube  retroreflector  design  speci¬ 
fies  a  tolerance  on  the  return-beam  direction  (up  to  3  arcsec  ss  15  //rad)  rather  than  the 
dihedral  angle  tolerances.  The  dihedral  angle  tolerance  can  nevertheless  be  inferred  from 
the  reflected  angle  tolerance,  and  this  has  been  done  in  Appendix  B.  It  is  found  that  the 
corresponding  dihedral  angle  tolerance  is  ±1.2  arcsec  (or  ±5.8  //rad),  for  a  cornercube 
of  this  material  and  geometry.  Note  that  the  unit  of  arcsec  is  used  throughout  this  work 
for  specifying  the  size  of  small  angles,  as  is  usual  in  the  field  of  SLR.  The  approximate 
conversion  to  //rad  is  possible  simply  by  multiplying  by  a  factor  of  4.85. 

The  return  SLR  laser  beam  can  be  broadened  (or  narrowed)  in  different  directions  by 
changing  the  dihedral  angles.  Broadening  the  return  beam  will  lower  its  radiance,  but 
often  this  penalty  is  outweighed  by  the  benefit  of  mitigating  the  beam  deflection  away 
from  the  SLR  receiver  due  to  the  velocity  aberration,  discussed  below.  Often,  satellite 
retroreflectors  will  be  intentionally  adjusted  (known  as  spoiling)  in  this  way.  Typically, 
the  amount  of  spoiling  will  be  less  than  approximately  2  arcsec  [5].  The  spoiling  might 
be  applied  equally  to  all  the  dihedral  angles  for  a  more-or-less  rotationally  symmetri¬ 
cal  change  in  beam  width,  or  different  amounts  applied  to  each  angle  for  preferential 
broadening  in  some  directions  but  not  others.  Otsubo  [6]  gives  examples  of  spoiling  of 
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0.35  arcsec  that  broadens  (or  reduces,  depending  on  whether  it  is  applied  to  all  angles  or 
not)  the  beam  width  by  a  factor  of  two. 

The  retroreflectors  of  the  current  design  therefore  have  a  dihedral  angle  tolerance  that 
could  lead  to  spoiling  that  profoundly  affects  their  performance  (one  way  or  another)  in 
SLR.  In  the  present  work,  a  perfect  geometry  is  assumed  in  order  to  determine  a  base¬ 
line  attitude  control  strategy.  It  is  nevertheless  recommended  that  the  far-field  diffrac¬ 
tion  pattern  of  each  retroreflector  be  characterised  experimentally  in  order  to  determine 
its  optimal  attitude  control  strategy  due  to  the  potential  for  significant  differences  due 
to  the  actual  detailed  cornercube  geometry  of  the  retroreflectors.  An  industry-standard 
method  for  testing  space  retroreflectors  is  described  by  Boni  et  al.  [7]  and  by  Dell'Agnello 
et  al.  [8].  The  method  uses  a  laser-beam  profiler,  beam  splitter,  and  an  arrangement  of 
polarising  prisms  and  waveplates  to  measure  the  diffraction  pattern  resulting  from  uni¬ 
form  laser  beams  in  different  polarisation  states.  A  different  wavelength  laser  light  is 
used  in  those  works  from  the  present  SLR  station,  although  this  need  not  be  the  case.  If 
a  different  wavelength  is  used,  the  diffraction  pattern  will  be  different,  and  adjustments 
to  the  measured  reflectance  pattern  would  need  to  be  made  in  order  to  predict  the  actual 
SLR  reflectance  pattern.  This  could  be  done  using  the  simulation  methods  developed,  for 
example,  by  Arnold  [9]. 


4  Radar  link  equation 


The  number,  np,  of  photons  detected  at  an  optical  receiver,  from  a  laser  source  after  re¬ 
flection  from  a  (distant)  satellite,  is  often  calculated  using  an  expression  from  the  theory 
of  radar  called  the  radar  link  equation.  This  expression  includes  the  effects  of  losses  within 
the  transmitting  and  receiving  optics,  losses  during  the  journey  of  the  beam  through  the 
atmosphere  in  both  directions,  and  the  effects  of  beam  broadening  on  both  the  outward 
and  return  paths.  The  present  work  follows  the  methods  outlined  in  the  review  article  by 
Degnan  [5]. 

The  radar  link  equation  for  the  number  of  detected  photons  can  be  written  as 


np  = 


Et\  Gt  n 
he  ^  And2  And2 


ArrirhqTa  U  • 


(1) 


Here  the  first  term  accounts  for  the  number  of  photos  emitted  during  a  pulse  of  the  laser 
source  of  energy  Et  and  wavelength  A,  with  h  the  Plank  constant  and  c  the  speed  of  light 
in  a  vacuum.  The  efficiency,  rjt,  of  the  transmitting  optics  accounts  for  losses  within  the 
broadcasting  telescope. 

The  reduction  in  the  irradiance,  due  to  broadening,  of  the  beam  on  its  path  of  length  d 
(also  known  as  the  range)  to  the  retroreflector  is  given  by  the  term  Gt  /  {^nd?)  where  the 
gain,  Gt,  is  a  function  of  beam  divergence  and  direction,  as  discussed  below.  Similarly, 
the  effect  of  the  divergence  of  the  light  reflected  from  the  satellite  (retroreflector)  is  de¬ 
scribed  by  the  subsequent  term,  where  a  is  the  optical  cross  section  of  the  retroreflector. 
The  optical  cross  section  is  due  to  the  finite  retroreflector  area  presented  to  the  incoming 
beam  and  the  effects  of  diffraction  of  the  light,  and  so  depends  on  both  the  orientation  of 
the  retroreflector,  and  the  direction  of  the  receiver. 
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The  return  path  to  the  receiver  is  assumed  to  be  of  equal  length,  d,  to  that  of  the  outward 
path.  During  the  double  crossing  of  the  Earth's  atmosphere,  a  proportion  of  the  light  will 
be  lost  through  absorption  and  scattering.  This  is  accounted  for  by  the  transmittances  Ta 
and  Tc  through  the  atmosphere  and  cirrus  clouds  respectively. 

The  receiver  will  have  receiving  optics  of  area  Ar  and  receiving  optics  efficiency  rjr,  while 
the  efficiency  of  the  detector  itself  is  r/q. 

Expressions  for  the  gain,  Gt,  and  the  retroreflector  optical  cross  section,  a,  were  given  in 
S&P  [1],  and  are  summarised  below.  Values  for  the  atmospheric  transmittance,  Ta,  and 
the  cirrus-cloud  transmittance,  Tc,  were  also  calculated  for  a  range  of  likely  scenarios. 
The  present  work  largely  concerns  the  retroreflector  optical  cross  section,  or  gain,  a,  its 
variation  with  the  retroreflector  orientation,  and  the  effect  this  has  on  the  final  photon 
count  at  the  receiver. 


4.1  Transmitter  gain 


The  transmitter  gain  depends  on  the  propagation  characteristics  of  the  laser  beam;  which 
in  turn  depend  on  the  transmitter  optics,  aperture  size,  beam  divergence  and  beam  pro¬ 
file.  These  are  considered  below  for  the  present  application  of  satellite  laser  ranging 
where  the  laser  beam  has  a  top-hat  beam  profile.  In  this  case  the  Fresnel-Kirchhoff  diffrac¬ 
tion  integral  can  be  used  to  derive  an  expression  for  the  gain  [10, 11].  If  the  aperture  inner 
and  outer  radii  are  respectively  b  and  a,  and  if  the  radius  of  curvature  of  the  wavefront  at 
the  aperture  is  again  R,  then  the  gain 


Gt 


exp[i/3it]  Jq(X\/u)  du 


(2) 


where  Jo  is  a  Bessel  function  of  the  first  kind,  u  is  the  normalised  radial  distance  on  the 
laser  aperture,  7  =  b/a,  X  =  ka  sin  0\,  Q\  is  the  angle  from  the  central  axis  of  the  beam — in 
the  present  case,  the  pointing  error  of  the  laser — and 


P 


ka 2  / 1  1  \ 

~Y  \d  +  Rj  ' 


(3) 


The  (slant)  range,  d,  is  calculated  assuming  a  circular  satellite  orbit  of  altitude  hs  = 
600  km,  about  the  Earth  of  radius  Re  =  6.378  x  103  km.  For  a  pointing  angle  of  0  from 
the  zenith,  the  range  is 

d  =  \J R%  cos2  4>  +  2 Rehs  +  h?s  —  Re  cos  </>  (4) 


from  simple  trigonometry.  Note  that  no  account  of  atmospheric  refraction  has  been  made 
in  this  approximation.  Typically,  the  error  is  not  greater  than  10  m  [5]  and  so  can  be  safely 
neglected  when  calculating  return  pulse  strength. 

If  the  beam  divergence  is  again  denoted  9d,  then  the  wavefront  radius  of  curvature  at  the 
aperture  can  be  calculated  using  straightforward  trigonometry  as 


R  = 


a 

tan  y 


(5) 
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Figure  3:  Transmitter  gain,  Gt{9 i)  with  angle,  0  \ ,  of  the  laser  beam  direction  error.  Values 
of  the  beam  pointing  accuracy  metrics  (the  relative  and  absolute  rms  pointing  accuracies) 
are  shown  as  vertical  dashed  grey  lines.  The  dashed-dotted  red  line  shows  the  angle  of 
stationary  gain,  6\  =  0.085  arcsec.  The  gain  for  a  selection  from  the  allowed  range  of 
satellite  zenith  angles,  <j>,  and  therefore  distances,  d,  are  shown. 


A  positive  value  of  R  corresponds  to  a  diverging  beam,  while  a  negative  value  corre¬ 
sponds  to  a  converging  beam. 


The  resulting  gain  is  plotted  as  a  function  of  beam  direction  in  Figure  3  for  a  range  of 
pointing  angles  (and  therefore  of  satellite  ranges).  The  relative  and  absolute  rms  pointing 
accuracies  of  the  station  are  0.1  arcsec  and  2  arcsec,  respectively,  are  indicated  in  the 
figure.  The  magnitude  of  the  central  bright  spot  (9 1  <  0.08  arcsec)  and  first  lobe  (6\  < 
0.2  arcsec)  can  be  seen  to  vary  by  a  factor  of  approximately  10  and  2,  respectively,  over 
the  range  of  target  distances.  The  rest  of  the  distribution  is  nevertheless  relatively  stable. 
The  effects  of  atmospheric  turbulence  on  coherent-beam  propagation  is  too  complex  to  be 
considered  in  any  detail  in  the  present  work.  Degnan  [5]  states  that  the  resulting  beam¬ 
broadening  is  likely  to  be  minimal  in  most  cases,  being  not  greater  than  a  few  percent. 
Beam  wander  and  scintillation  of  the  beam  will  be  apparent  as  fluctuations  in  the  signal 
strength,  without  a  change  in  the  mean  value.  The  values  of  the  gain  at  a  beam  pointing 
error,  9\,  of  0.085  arcsec  can  be  seen  not  to  vary  with  satellite  direction,  &,  from  the  zenith. 
This  value  of  beam  pointing  error  is  close  to  the  relative  rms  pointing  accuracy.  So  the 
transmitter  gain  at  this  angle  (GY  (0.085)  =  1.25  x  1010)  is  taken  as  a  representative  worst- 
case  value  in  the  remainder  of  the  present  work. 
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4.2  Retroreflector  optical  cross-section 


As  shown  in  S&P  [1],  the  retroreflector  optical  cross-section,  a,  in  the  far-field  limit, 
rcc/  d  — >  0,  can  be  written 


(7  = 


47tA2c  2,/i  (krcc  sin  9r) 
^  A2  krcc  sin  9r 


(6) 


when  the  retroreflector  is  oriented  normal  to  the  incident  light,  and  where  6r  is  the  re¬ 
ceiver  orientation  error  (see  Figure  1),  p  is  the  retroreflector  reflectance  (or  efficiency),  J\ 
is  a  Bessel  function  of  the  first  kind,  k  is  the  propagation  number  2tt/A,  and  Acc  =  irrf.c  is 
the  area  of  the  retroreflector  circular  aperture  of  radius  rcc.  The  optical  cross-section  for 
on-axis  (that  is,  when  9r  =  0)  reflected  light  will  therefore  be 


a 


er= o 


4vr 

P^2A 


2 

cc 


(7) 


which  can  be  seen  to  be  proportional  to  the  square  of  the  retroreflector  aperture  area. 
Light  at  other  reflection  angles  is  distributed  according  to  the  Airy  diffraction  pattern  as 
described  above  in  equation  (6). 

For  on-axis  reflected  light,  when  the  retroreflector  is  a  corner  cube,  and  is  not  oriented 
normal  to  the  incident  light  but  at  an  angle  9t  from  the  normal,  the  area  Acc  in  the  above 
expression  (7)  must  be  reduced  by  the  factor  [12,  5]  (and  see  Appendix  A) 


K(9i) 


2 

7T 


[sin  1  p  —  (l / rcc)p  tan  \9'\]  cos  9i, 


(8) 


where  9’  is  the  refracted  angle  in  the  corner  cube  given  by  Snel's  law  of  refraction 


rising  =  sin#*, 


(9) 


where  n  is  the  cube  index  of  refraction,  and  where 


p  =  \J  1  -  (l /fee)2  tan2  9'  (10) 

is  the  normalised  reduced  aperture  radius  along  the  tilt  direction.  The  variation  of  the 
optical  cross-section  (on  axis)  with  retroreflector  angle  is  therefore  k2  ( 9,, ) .  This  function  is 
shown  plotted  in  Figure  4  for  the  present  case  of  a  retroreflector  made  of  BK7  crown  glass. 
Its  values  are  tabulated  in  Table  1.  The  cornercube  ratio  l/rcc  is  an  important  parameter 
of  the  retroreflector.  It  is  the  ratio  of  the  depth,  l,  of  the  cornercube  to  its  aperture  radius, 
rcc.  For  a  fully  circular  aperture,  the  ratio  must  be  at  least  y/2,  but  is  often  larger  than  this 
in  order  to  avoid  thin  edges.  The  apparent  aperture  area  decreases  with  tilt  angle  most 
slowly  in  cornercubes  with  this  minimum-value  ratio 

(Z/rcc)min  -  y/2  (11) 

as  shown  in  Figure  4.  When  this  ratio  is  larger,  as  in  the  present  case,  the  fall-off  is  more 
rapid,  and  such  retroreflectors  are  more  sensitive  to  the  effects  of  tilting. 

More  generally  the  optical  cross-section  off-axis  is  broadened  in  one  direction  by  the  nar¬ 
rower  effective  aperture  in  the  direction  of  the  tilt.  In  this  case  the  optical  cross-section 
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Figure  4:  On-axis  retroreflector  optical  cross-section  variation,  k2(0,)  with  angle,  0U  of  the 
retroreflector  axis  from  the  incident  beam  direction.  The  plot  shows  the  case  of  a  solid 
retroreflector  constructed  of  BK7  crown  glass  with  refractive  index  n  =  1.519.  Values  of 
the  angle  6^  are  shown  where  the  gain  has  reduced  to  80%,  50%  and  10%  of  its  maximum 
value. 


is  more  difficult  to  calculate  because  of  the  more  complex  aperture  shape  in  the  diffrac¬ 
tion  integral.  The  method  of  calculation  is  described  in  Appendix  A  and  the  resulting 
expression  for  the  optical  cross-section  is 

An t  A 2 

v  =  P^^W2(x,y)  (12) 

Az  7TZ 

where  W (x,  y)  is  defined  by  equation  (All)  but  calculated  using  equations  (A16)  and  (A17) 
On  axis,  note  that 

W(0, 0)  =  7i7«(0i)  (13) 

as  expected.  Off  axis,  the  effects  of  the  tilted  retroreflector  are  illustrated  in  Figure  5. 
The  optical  cross-section  can  be  seen  to  broaden  in  the  same  direction  that  the  effective 
aperture  width  is  reduced,  as  discussed  in  Appendix  A.  The  overall  magnitude  of  the 
optical  cross-section  is  reduced,  in  accordance  with  the  variation  of  n2(9i),  in  the  on-axis 
(9r  =  0)  direction. 

4.2.1  Velocity  aberration 

As  the  receiver  is  co-incident  with  the  transmitter  optics,  one  might  expect  the  appro¬ 
priate  value  for  the  reflected  light  angle,  9r,  in  the  radar  equation  to  be  zero.  Indeed, 
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Table  1:  Retroreflector  optical  cross-section  variation,  K2{9i),  tabulated  as  a  function  of 
angle,  9i,  of  the  retroreflector  axis  from  the  incident  beam  direction.  In  this  case  the  retrore¬ 
flector  is  made  of  solid  BK7  crown  glass  with  refractive  index  n  =  1.519. 


0i(°) 

K2 

Av2 

W 

K2 

W 

K2 

W 

K2 

0.0 

1.000 

6.0 

0.695 

12.0 

0.441 

18.0 

0.248 

24.0 

0.118 

0.5 

0.973 

6.5 

0.672 

12.5 

0.422 

18.5 

0.235 

24.5 

0.110 

1.0 

0.946 

7.0 

0.649 

13.0 

0.404 

19.0 

0.222 

25.0 

0.102 

1.5 

0.920 

7.5 

0.626 

13.5 

0.387 

19.5 

0.210 

25.5 

0.095 

2.0 

0.894 

8.0 

0.604 

14.0 

0.369 

20.0 

0.198 

26.0 

0.088 

2.5 

0.868 

8.5 

0.582 

14.5 

0.353 

20.5 

0.187 

26.5 

0.081 

3.0 

0.842 

9.0 

0.561 

15.0 

0.336 

21.0 

0.176 

27.0 

0.075 

3.5 

0.817 

9.5 

0.540 

15.5 

0.321 

21.5 

0.165 

27.5 

0.069 

4.0 

0.792 

10.0 

0.519 

16.0 

0.305 

22.0 

0.155 

28.0 

0.063 

4.5 

0.767 

10.5 

0.499 

16.5 

0.290 

22.5 

0.145 

28.5 

0.058 

5.0 

0.743 

11.0 

0.479 

17.0 

0.276 

23.0 

0.136 

29.0 

0.053 

5.5 

0.719 

11.5 

0.460 

17.5 

0.262 

23.5 

0.127 

29.5 

0.048 

this  would  be  correct  if  the  target  (satellite)  were  stationary  relative  to  the  transmit¬ 
ter/  receiver  station.  Yet  the  relative  velocity  of  the  satellite  is  in  general  non-zero,  and 
this  causes  a  small,  but  in  many  cases  significant,  change  in  the  apparent  direction  of  the 
laser  beam — known  as  the  velocity  aberration  [5, 13].  Denoted  by  a,  the  aberration  causes 
a  change  in  the  receiver  angle  0r  =  a  as  illustrated  in  Figures  6  and  1.  As  shown  below, 
the  magnitude  of  the  aberration  can  be  as  large  as  approximately  11  arcsec;  comparable 
to  the  width  of  the  reflected  beam,  so  the  effect  is  significant. 

The  velocity  of  the  satellite  is  sufficiently  large  (ss  7.5  km/s)  to  affect  appreciably  the 
apparent  angle  at  which  the  laser  beam  strikes  it.  If  the  component  of  the  satellite  rela¬ 
tive  velocity,  v,  that  is  orthogonal  to  the  direction  to  the  station  is  denoted  v1,  then  the 
change  in  angle  is  [14]  fair  1  (t;//c)  by  the  classical  (Bradley)  approximation,  which  holds 
when  v'  -C  c  where  c  is  the  speed  of  light.  The  reflected  beam  will  undergo  a  similar 
perceived  rotation  in  the  station  frame  of  reference,  and  so  the  total  angular  shift — or 
velocity  aberration — will  be 


a  =  2  tan  1 


2v' 

c 


(14) 


This  first-order  approximation  is  sufficiently  accurate  for  satellite  speeds,  where  the  second- 
order  terms  have  been  shown  [13]  to  be  approximately  10~4  times  smaller. 

The  speed  of  a  satellite  in  a  circular  orbit  can  be  simply  calculated  from  its  altitude,  hs,  as 


v  = 


1  Rig 

Re  +  hs 


(15) 


where  Re  is  the  radius  of  the  Earth  and  g  is  the  acceleration  due  to  gravity.  It  can  be 
shown  [15,  5]  that  the  velocity  aberration  is 

a(4>,  uj)  =  am\J  cos2  oj  +  T2(</>)  sin2  u  (16) 
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Diffraction  patterns  for  tilted  retroreflector,  rcc  =  6.350  mm 
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Figure  5:  Normalised  optical  cross-section,  cr/tTo,  as  a  function  of  reflected  angles  9r’v^ 
defined  by  ( x ,  y)  =  <1  tan  ()l:,x'v) .  False-colour  plots  corresponding  to  each  of  four  tilt  angles, 
9i,  of  the  retroreflector  about  an  axis  parallel  to  the  y-axis  are  shown.  Broadening  of  the 
diffraction  pattern  along  the  cc-axis  can  be  seen,  accompanied  by  an  overall  reduction  in 
magnitude  according  to  k2(A)  on  axis  (9lf  '1  =  9 =  0). 


where  (j)  is  the  transmitter  pointing  angle  from  the  zenith,  the  quantity  F  varies  with 
pointing  direction  and  is  given  by 


m 


/  Re  sin  <f>\2 

V  Re  +  hs  ) 


(17) 


and  u)  is  the  azimuthal  angle  of  the  satellite  velocity  in  the  plane  orthogonal  to  the  posi¬ 
tion  vector,  r,  of  the  satellite  relative  to  the  station.  The  angle  uj  can  be  shown  to  be 

ui  =  cos_1[(f  x  s)  ■  b]  (18) 


where  s  is  the  position  vector  of  the  satellite  relative  to  the  centre  of  the  Earth,  and 
indicates  a  unit  vector.  The  velocity  aberration  therefore  varies  between  the  two  extremes 

Chnax  —  Cbn  (19) 
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z 


Figure  6:  Diagram  illustrating  the  beam  directions  in  both  the  ground-  and  satellite-based 
frames  of  reference  that  account  for  the  velocity  aberration,  a,  in  the  reflected  light  distri¬ 
bution.  Figure  taken  from  [13]. 


and 


Chnin  — 


(20) 


for  a  given  pointing  angle  <f>,  and  where  the  maximum  aberration,  am,  depends  on  the 
speed  and  can  be  expressed  as 

2v 

am  =  — •  (21) 

c 

If  the  rotation  of  the  Earth — and  hence  the  velocity  of  the  station,  v — is  taken  in  to  ac¬ 
count,  then  for  a  general  orbit  path  the  two  extremes  become 


and 


2{v  +  v)  2v 

Chnax  =  =  (%m  H  j 

c  c 


_  2(vT(</>)-v)  _  ~  2v 

^min  —  —  (0) 

C  C 


The  station  speed  depends  on  its  latitude,  <p,  and  is,  from  simple  geometry. 


(22) 

(23) 


27t Re  cos  (p 
Tday 


(24) 


where  Tday  is  the  length  of  one  day 

These  limits  to  the  range  of  the  velocity  aberration  for  a  satellite  in  a  600  km  orbit,  and  a 
station  at  Mount  Stromlo,  are  plotted  in  Figure  7  as  a  function  of  pointing  angle.  The  max¬ 
imum  velocity  aberration,  am,  is  the  same  at  all  pointing  angles  and  is  equal  to  10.92  arc- 
sec  (52.9  //rad).  The  minimum  varies  with  pointing  angle,  from  its  largest  at  zenith  down 
to  amin  =  4.4  arcsec  (21.3  //rad)  at  75°  from  the  zenith. 

This  range  of  velocity  aberrations  applied  to  the  current-design  retroreflector  optical  cross 
section  can  be  seen  in  Figure  8.  The  unity  photon-count  threshold  is  also  shown,  and  was 
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Figure  7:  Velocity  aberration,  a,  range  (shaded)  as  a  function  of  pointing  angle  for  an  orbit 
height  of  600  km.  The  range  includes  the  effects  due  to  a  circular  satellite  orbit,  and  the 
speed  of  the  station  due  to  the  rotation  of  the  Earth. 


calculated  by  inverting  the  radar  link  equation  (1)  with  the  photon  count  np  =  1,  and 
worst  expected-case  parameter  values  used  for  the  transmitter  gain,  transmittance  and 
range;  that  is,  9\  =  0.1  arcsec  and  (!)  =  75°.  This  threshold  will  therefore  be  lower  in  more 
favourable  circumstances.  It  can  be  seen  from  this  figure  that  the  presence  of  the  larger 
velocity  aberrations  will  cause  the  minimum  of  the  diffraction  pattern  of  the  retroreflector 
directed  at  the  station  to  coincide  with  the  receiver,  and  so  the  photon  count  to  fall  below 
the  required  unity  threshold. 

There  is  more  than  one  way  to  address  this  problem.  Very  often  retroreflectors  are  spoiled 
to  broaden  the  reflected  light  preferentially  to  the  angles  9r  close  to  the  range  of  expected 
velocity  aberration  a.  Spoiling  is  the  term  used  to  describe  a  cornercube  retroreflector 
design  with  faces  meeting  at  angles  slightly  perturbed  (<  2  arcsec)  from  the  usual  right- 
angles  [5,  6].  This  approach  is  not  possible  in  the  present  project. 

An  alternative  is  to  use  a  retroreflector  with  a  smaller  aperture,  and  hence  a  broader 
diffraction  pattern  whose  minima  fall  outside  the  range  of  velocity  aberration  angles.  The 
diffraction  pattern  of  just  such  a  retroreflector  is  plotted  in  Figure  8.  Although  this  does 
increase  the  optical  cross  section  above  the  critical  threshold  where  required,  it  reduces 
its  value  considerably  in  other  regions  and  places  greater  constraints  on  the  range  of 
acceptable  tilt  angles;  requiring  greater  satellite  attitude  control. 

The  final,  and  preferred,  option  is  to  broaden  the  diffraction  pattern  of  the  existing  retrore¬ 
flector  by  tilting  it  to  achieve  a  narrower  effective  aperture  in  the  same  direction  as  the 
velocity  aberration.  If  the  normal  to  the  retroreflector  aperture  is  tilted  as  an  angle  0, 
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Figure  8:  Normalised  retroreflector  optical  cross  section,  <t/<to,  as  a  function  of  receiver 
angle,  9r,  for  a  range  of  retroreflector  tilt  angles  and  two  aperture  sizes.  The  range  of 
applicable  receiver  angles  due  to  the  velocity  aberration  is  shaded  grey.  The  unity  photon- 
count  threshold  is  shown  shaded  yellow. 

away  from  the  station,  and  in  the  direction  of  the  projected  relative  satellite  velocity 


(25) 


then  the  diffraction  pattern  is  broadened  in  the  direction  of  the  velocity  aberration.  In 
this  way  the  optical  cross  section  in  increased  above  the  threshold  where  (and  when) 
required. 


The  optimal  tilt  angle,  9f  \  can  be  found  for  any  given  velocity  aberration  by  numeri¬ 


cally  maximising  the  optical  cross  section  with  respect  to  0*.  The  optimal  angle  will,  in 


general,  vary  with  the  velocity  aberration,  a.  The  resulting  function,  6^°\a),  is  plotted 


in  Figure  9.  For  velocity  aberrations  a  below  approximately  5.5  arcsec,  any  increase  in 
the  optical  cross  section  due  to  beam  broadening  with  tilt  is  insufficient  to  compensate 
for  the  accompanying  reduction  in  the  optical  cross  section  due  to  the  reduced  effective 
retroreflector  area  presented  to  the  laser  beam.  In  these  cases,  directing  the  retroreflector 
aperture  directly  at  the  station  is  optimal;  that  is,  when  Oi  =  0.  By  contrast,  at  higher 
velocity  aberration  values,  tilting  the  retroreflector  is  advantageous.  For  example,  when 
the  velocity  aberration  is  its  highest  (for  the  present  orbit)  a  =  amax  =  10.92  arcsec,  the 
optimal  tilt  angle  (9?-  =  20.5°,  and  the  range  of  tilt  angles  for  which  the  optical  cross 

section  exceeds  the  unity  photon-count  threshold  is  quite  wide;  namely,  5°  <  6i  <  39°. 
If  necessary,  a  tilt  angle  within  this  range  centred  at  0t  =  20.5°  could  be  maintained  for 
any  velocity  aberration  value,  and  the  resulting  optical  cross  section  would  remain  above 
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Figure  9:  Contour  plot  of  the  (normalised)  optical  cross  section  of  the  retroreflector  design 
of  the  present  project,  at  receiver  angle  8r  =  a,  the  velocity  aberration.  The  variation  of 
the  optical  cross  section  with  tilt  angle  in  the  direction  required  to  broaden  the  reflected 
beam  toward  the  receiver  is  shown.  The  red  curve  is  the  loci  of  optimal  tilt  angles  d[°’ (a) 
that  maximise  the  optical  cross  section  of  the  retroreflector.  Contours  of  the  optical  cross 
section  corresponding  to  ten  times  the  unity  photon-count  threshold  are  drawn  to  indicate 
the  range  of  acceptable  tilt  angles. 


the  threshold,  as  shown  in  Figure  8.  Yet  this  would  not  be  optimal  for  many  situations. 
If  possible,  tilting  by  the  optimal  angles  shown  by  the  curve  in  Figure  9  would  give  a 
higher  photon  count.  To  this  end,  the  calculated  values  of  the  optimal  tilt  angles  f)il°> (a) 
are  listed  in  Table  2. 


The  tolerance  for  the  direction  of  the  retroreflector  tilt  axis  is  quite  large,  as  shown  in 
Figure  10.  Denoting  the  azimuthal  error  of  the  tilt  axis  direction  by  ipt,  it  can  be  seen 
from  the  figure  that  if  the  tilt  angle  magnitude  is  optimal  (that  is,  0r  =  6^°\a))  then 
the  unity  photon-count  threshold  is  exceeded  for  all  expected  velocity  aberration  values 
when  ipt  <  70°.  This  wide  tolerance  for  the  unity  threshold  is  much  larger  than  the 
tolerance  required  to  meet  the  10  x  threshold.  At  the  higher  threshold  the  tolerance  on 
the  orientation  error  is  < pt  <  25°. 
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Table  2:  Values  of  the  optimum  retroreflector  tilt  angles,  9^°\  (in  degrees)  as  a  function  of 
velocity  aberration,  a,  (in  arcsec)  for  the  current  retroreflector  design. 


a 

of 

a 

of 

a 

of 

a 

of 

0.000 

0.0000 

7.000 

8.8510 

8.714 

15.030 

10.429 

19.280 

5.429 

0.0000 

7.143 

9.4735 

8.857 

15.442 

10.571 

19.584 

5.571 

0.9243 

7.286 

10.072 

9.000 

15.842 

10.714 

19.875 

5.714 

1.8964 

7.429 

10.649 

9.143 

16.229 

10.857 

20.161 

5.857 

2.8187 

7.571 

11.204 

9.286 

16.605 

11.000 

20.444 

6.000 

3.6959 

7.714 

11.740 

9.429 

16.970 

11.143 

20.703 

6.143 

4.5322 

7.857 

12.256 

9.571 

17.335 

11.286 

20.967 

6.286 

5.3304 

8.000 

12.756 

9.714 

17.682 

11.429 

21.235 

6.429 

6.0945 

8.143 

13.242 

9.857 

18.022 

11.571 

21.488 

6.571 

6.8273 

8.286 

13.715 

10.000 

18.351 

11.714 

21.722 

6.714 

7.5319 

8.429 

14.166 

10.143 

18.671 

11.857 

21.971 

6.857 

8.2109 

8.571 

14.605 

10.286 

18.983 

12.000 

22.199 

4.2.2  Point  ahead  ranging 


Pointing  the  transmitter  ahead  of  the  apparent  position  of  the  satellite  to  allow  for  the 
time  of  flight  of  the  laser  pulse  is  often  done  when  the  return  signal  needs  to  be  max¬ 
imised;  as  in  the  case  of  uncooperative  targets,  for  example.  The  size  of  the  point-ahead 
angle,  6\,  can  be  calculated  from  the  component,  v' ,  of  the  satellite's  relative  velocity  that 
is  orthogonal  to  the  position  vector,  r,  of  the  satellite  relative  to  the  station.  The  satellite 
relative  velocity,  v,  is  the  time  derivative  of  r,  and  (from  equation  (14)) 


v'  =  \v  —  (v  ■  r)r|  = 


ca 


(26) 


where  r  =  r/jr|,  and  a  is  the  velocity  aberration.  The  one-way  time  of  flight  will  be 


and  so  the  point-ahead  angle  can  be  written 


(27) 


a 
2  ’ 


(28) 


that  is,  half  the  velocity  aberration  angle. 

By  the  time  the  pulse  has  returned,  the  receiver  will  have  rotated  a  further  0\,  and  so  will 
be  at  an  angle  O2  =  a  to  the  returning  light  direction.  The  resulting  relative  change  in 
receiver  aperture  area  will  be  (1  —  cos  62)  ~  10-9,  which  is  negligible. 

These  calculations  neglect  the  effects  of  atmospheric  propagation.  Although  the  prin¬ 
ciples  of  the  above  discussion  still  hold,  in  practice  the  effects  of  atmospheric  refraction 
must  be  included  for  accurate  SLR.  These  calculations  are  beyond  the  scope  of  the  present 
work.  Further  details  of  a  standard,  and  long-standing,  shell  model  of  the  atmosphere 
known  as  the  Marini-Murray  model  are  given  in  the  review  article  by  Degnan  [5].  It  is 
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er  =  a,  e{  =  6t(o),  rcc  =  6.350  mm,  (t/rj2  =  3.41 


a  (arcsec) 


10 


10 


10 


Figure  10:  Contour  plot  of  the  (normalised)  optical  cross  section  of  the  retroreflector  design 
of  the  present  work,  at  receiver  angle  9r  =  a,  the  velocity  aberration.  The  retroreflector 
tilt  angle  e,  =  the  optimum  value  for  each  value  of  the  velocity  aberration,  but 

the  direction  of  the  tilt  axis  is  in  error  by  the  azimuthal  angle  The  unity  photon-count 
threshold  contour  is  also  plotted. 


simply  noted  in  passing  that  this  particular  model  requires  values  of  physical  measures 
of  the  atmosphere — such  as  surface  relative  humidity,  temperature,  and  pressure — to  im¬ 
plement. 

It  should  be  noted  that  the  point-ahead  principle  does  not  correct  the  velocity  aberration, 
despite  the  direct  relationship  to  it,  since  the  direction  of  the  part  of  the  beam  incident 
on  the  retroreflector  is  unchanged,  and  hence  the  reflected  beam  is  still  returned  at  the 
aberrated  angle  a.  Nevertheless,  it  does  ameliorate  some  of  its  effect  by  maximising  the 
radiance  of  the  beam  striking  the  retroreflector.  A  point-ahead  strategy  is  assumed  in  the 
photon-budget  calculations  in  the  present  work. 


4.3  Atmospheric  and  cirrus-cloud  transmittance 

The  single-transit  atmospheric  transmittance,  Ta,  values  were  calculated  using  the  MOD- 
TRAN  5  [16]  computer  software.  The  built-in  mid-latitude  summer  and  winter  atmo¬ 
spheric  models  were  used.  The  visibility  component  was  set  to  the  Rural  Aerosol  model 
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Figure  11:  One-way  atmospheric  transmittance,  Ta,  calculated  using  MODTRAN  5  [16] 
for  light  of  wavelength  532  run  and  a  bandwidth  of  1  run.  Calculations  were  for  the  mid¬ 
latitude  summer  and  winter  atmospheric  models,  rural  aerosol  visibility  of  23  km,  and  a 
ground  altitude  of  770  m  above  sea  level. 


with  visibility  of  23  km.  The  ground  altitude  above  mean  sea  level  was  chosen  to  be 
770  m,  matching  the  Mount  Stromlo  station  altitude.  Light  of  wavelength  532  nm  and 
with  a  bandwidth  of  1  nm  was  modelled. 

The  resulting  transmittance  values  for  the  useful  range  of  zenith  pointing  angles  are  plot¬ 
ted  in  Figure  11.  At  any  pointing  angle,  the  summer  atmosphere  model  gave  a  slightly 
smaller  transmittance  than  the  winter  model  at  the  same  angle.  The  summer  model  trans¬ 
mittance  values  were  therefore  used  in  establishing  the  retroreflector  tilt-angle  tolerances. 

The  other  transmittance  parameter  is  the  single-transit  cirrus-cloud  transmittance,  Tc. 
Degnan  [5]  reported  that  this  could  be  well  represented  by  a  value  of  approximately 
0.8.  This  value  is  the  mean  cirrus-cloud  transmittance  at  zenith  when  cirrus  clouds  are 
present,  and  has  been  found  to  be  exceeded  75%  of  the  time.  The  mean  transmittance  is 
lower  at  other  angles;  as  low  as  0.1  at  70°  pointing  angle.  Yet  it  is  found  that  no  cirrus 
cloud  is  present  50%  of  the  time,  corresponding  to  a  transmittance,  Tc  =  1.  Accordingly, 
the  representative  value  of  Tc  =  0.8  is  used  in  the  present  calculations. 


5  Photon  budget 

The  results  of  a  calculation  of  the  photon  count  for  the  current  system  based  on  estimates 
of  alignment  tolerances  are  shown  in  Table  3.  The  table  also  lists  the  relevant  parameter 
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values  used.  The  retroreflector  efficiency  value,  p,  is  based  on  the  assumption  of  three  re¬ 
flections  of  the  light  from  the  silvered  BK7  surfaces,  for  which  a  single-surface  reflectance 
of  0.95  is  assumed  [2],  The  other  parameter  values  were  discussed  above. 

The  first  two  cases  of  the  table  show  the  photon  counts,  np,  when  the  satellite  is  at  the 
zenith  (0  =  0),  and  its  velocity  aberration  is  at  its  smallest  for  this  pointing  angle;  that  is, 
a  (=  6r)  =  9.88  arcsec.  With  no  tilting  compensation,  the  photon  count  np  =  625.0  which 
is  adequate  for  pulse  detection.  The  photon  count  increases  approximately  25-fold,  which 
is  substantial,  on  applying  the  optimum  tilt  angle  of  9i  =  18.2°. 

Cases  3-5  show  the  photon  counts  in  the  more  marginal  situation  when  the  pointing  an¬ 
gle,  cf)  =  75°  from  the  zenith.  The  velocity  aberration  varies  between  4.37  arcsec  and 
10.92  arcsec,  and  without  any  tilt  compensation  the  photon  counts  are  respectively  184.1 
and  0.279 — the  latter  being  barely  sufficient  for  pulse  detection,  even  in  darkness.  Ap¬ 
plication  of  the  optimum  tilt  angle  9\  =  20.5°  in  the  latter  case  yields  a  photon  count  of 
21.22,  which  is  now  adequate  for  pulse  detection. 


6  Conclusion 


In  conclusion,  it  has  been  found  that  the  retroreflector  with  ideal  geometry  should  be 
tilted  away  from  the  SLR  station  in  order  to  mitigate  the  effects  of  the  velocity  aberration. 
The  optimum  tilt  angle  varies  with  the  size  of  the  velocity  aberration.  At  the  maximum 
expected  velocity  aberration  of  10.92  arcsec,  a  retroreflector  tilt  of  (20.5  ±  17)°  is  required 
to  exceed  the  unity  photon-count  threshold,  with  a  tilt  of  20.5°  being  optimal.  The  axis 
of  the  rotation  is  perpendicular  to  both  the  satellite  (relative)  velocity  and  the  direction 
to  the  SLR  station.  The  same  tilt  angle  (and  tolerance)  can  be  used  at  lower  velocity 
aberration  values  and  still  exceed  the  detection  threshold.  Nevertheless,  a  reduced  tilt 
angle  (according  to  the  values  in  Table  2  and  plotted  in  Figure  9)  will  give  optimal  pho¬ 
ton  counts.  The  most  restricted  range  of  acceptable  tilt  angles  occurs  when  the  velocity 
aberration  is  highest,  and  so  this  is  the  basis  for  the  attitude  control  tolerance.  Therefore, 
to  just  achieve  the  unity  photon-count  threshold  in  all  situations,  the  attitude  should  be 
within  17°  of  the  optimum  tilt  angle. 

The  retroreflector  manufacturer's  tolerance  on  the  reflected  beam  direction  of  within 
3  arcsec  has  been  found  to  be  equivalent  to  a  tolerance  of  ±1.2  arcsec  on  the  cornercube  di¬ 
hedral  angles.  This  amount  of  variation  encompasses  a  wide  range  of  far-field  diffraction 
patterns  that  deviate  significantly  from  that  of  the  ideal  geometry  considered  above.  It  is 
therefore  recommended  that  each  retroreflector  be  experimentally  characterised  to  deter¬ 
mine  its  diffraction  pattern  (and  its  optical  cross  section),  and  thence  its  associated  opti¬ 
mal  attitude-control  strategy  based  on  the  methods  described  in  the  present  work.  The 
characterisation  of  the  diffraction  pattern  could  be  a  direct  measurement  as  described,  for 
example,  by  Boni  et  al.  in 
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Table  3:  Parameter  values  used  in  calculating  the  detected  photon  number,  np.  Only  changes  in  parameter  values  have  been  recorded  in 
moving  between  adjacent  columns  from  left  to  right.  Values  marked  "1"  are  simple  estimates. 

Parameter  Case  1  Case  2  Case  3  Case  4  Case  5 
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Appendix  A  Tilted  cornercube  retroreflector  far-field 

diffraction  pattern 


The  optical  cross  section  of  a  retroreflector  illuminated  by  a  uniform,  coherent  light  across 
its  aperture  is 

a  =  Fnd2  p  U2(x,  y)  (Al) 

where  p  is  the  reflectance  of  the  retroreflector  and  U(x,  y)  is  the  electric  field  amplitude 
distribution  from  the  retroreflector  on  a  plane  with  coordinates  ( x ,  if),  at  the  observation 
distance  d.  In  the  present  case,  the  distance  d  is  much  larger  than  the  retroreflector  aper¬ 
ture  (a  circle  of  radius  rcc),  and  so  U(x,  y)  can  be  represented  accurately  by  the  far-field 
(Fraunhofer)  diffraction  distribution  [17] 


y)  =  exp  [id2  +  (^2  +  y2)] } 


exp 


k 


didr)  (A2) 


where  the  integral  is  over  the  effective  aperture,  f 1,  of  the  retroreflector,  with  Cartesian 
pupil-plane  coordinates  (£,  fj),  and  where  the  £-axis  is  parallel  to  the  x-axis.  The  wave¬ 
length  of  the  monochromatic  light  is  A,  and  k  =  ^  is  the  propagation  number. 

The  shape  of  the  aperture,  and  hence  the  limits  of  integration,  are  affected  by  the  tilt  of  the 
retroreflector.  In  the  case  of  a  circular  aperture  tilted  by  angle  Oi  from  normal  incidence 
of  the  light,  the  aperture  width  in  one  dimension — without  loss  of  generality  along  the 
£-axis — is  reduced  by  a  factor  of  cos#;.  Additionally,  in  cornercube  retroreflectors  there 
is  a  lateral  displacement  [9,  18]  of  the  image  of  the  aperture  which  in  turn  leads  to  a 
truncation  of  the  effective  aperture.  If  the  depth  of  the  cornercube  retroreflector  is  /,  and 
if  the  retroreflector  front  surface  is  not  recessed,  then  the  lateral  displacement  is  l  tan \0'\, 
where  O'  is  the  refracted  angle  of  the  light  ray  to  the  surface  normal  in  the  cornercube. 
This  angle  can  be  calculated  from  Snel's  law,  and  is 


O'  =  sin  1 


(A3) 


where  n  is  the  refractive  index  of  the  solid  cornercube  material.  In  the  case  of  a  hollow 
cornercube,  n  =  1  and  so  O'  =  0, .  Otherwise  O'  is  in  general  less  than  the  tilt  angle 
Oi.  The  resulting  expressions  describing  the  shape  of  the  effective  aperture  of  the  tilted 
cornercube  retroreflector  are  then 


f  =  ±A(??),  \r]\  <  p 

(A4) 

where 

F(fj)  =  cos  Oi  ^ \A'cC  —  r)2  —  /  tan \0'  j 

(A5) 

and 

p  =  F-1(0)  =rcc\j\-  {l/rcc)2  tan2  O'. 

(A6) 

Note  that 

in  the  untilted  case,  0,  =  O'  =  0  and 

F(v)  k=o  =  VrcC~v2 

(A7) 
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Figure  Al:  The  apparent  aperture  shapes  (solid  lines)  of  a  unit-radius,  solid  N-BK7  cor- 
nercube  retroreflector  at  a  number  of  different  tilt  angles,  0,.  The  retroreflector  has  been 
tilted  about  an  axis  parallel  to  the  y-axis.  The  outline  of  the  circular  aperture,  foreshort¬ 
ened  due  to  the  retroreflector  tilt  are  drawn  in  dashed  lines.  The  effective  apertures  are 
further  restricted  due  to  the  lateral  displacement  of  the  reflected  light,  by  an  amount  that 
increases  with  tilt  angle.  The  effective  aperture  is  equal  to  the  actual  aperture  only  in  the 
untilted  ( 9i  =  0°)  case,  shown  in  black. 


is  the  equation  of  a  circular  aperture  of  radius  rcc,  as  required.  Figure  Al  shows  plots 
of  the  effective  aperture  outlines  (£  =  ±F(fj))  of  a  unit-radius  retroreflector  at  various 
amounts  of  tilt.  Also  shown  are  the  corresponding  projected  aperture  outlines — the  con¬ 
siderable  effect  of  the  image  lateral  displacement  on  the  effective  aperture  can  be  clearly 
seen. 

The  evaluation  of  the  above  diffraction  integral  can  be  greatly  simplified  by  a  change  in 
variables.  Firstly,  let 

(x,y)  =  ^p(x,y)  (A8) 

and  then  make  use  of  the  normalised  pupil-plane  coordinates 

fa,v)  =  —fa,v).  (A9) 

Tec 

With  these  changes  to  the  variables,  the  above  expression  for  the  electric-filed  distribution 
at  the  observation  plane  may  be  written 

U2(x,y)  =  U(x,y)U*(x,y)  =  -j^W2(x,y)  (A10) 


UNCLASSIFIED 


23 


DST-Group-TR-3 1 72 


UNCLASSIFIED 


-0.6  -0.4  -0.2  (L0  0/2  04  oA 


e 

Figure  A2:  The  effective  aperture  (black)  of  a  unit-radius  aperture  tilted  9i  =  20° 
((I /rCc)2  =  3.41,  n  =  1.519),  and  an  example  set  of  N  =  8  equal-height  trapezoids  used 
to  approximate  the  effective  aperture  in  the  numerical  evaluation  of  the  diffraction  inte¬ 
gral.  In  practice,  a  much  larger  set  of  trapezoids  is  used. 


where 


and  where 


ni  rF(rj) 

W(x,y)=  /  exp  [— *  (£z  +  r/y)]  d£  drj 

J-p  9-F{r{) 


d  =  n! rcc  = 


F{rj)  =  F(y)/rcc  =  cos  Oi 


-  (l/rcc)2  tan2  6' 

(v/l  -V2  -  i  tan|0'|)  . 


Note  that  when  (x,  y)  =  (0, 0) — that  is,  the  receiver  orientation  error,  0r 
can  be  evaluated  analytically,  and  so  may  be  written 


W  (0, 0)  =  2  cos  6j, 


sin  1  n  —  J  y  t  an  1 6' 


=  Fn{0i). 


(All) 

(A12) 
(A13) 
0 — the  integrals 


(A14) 


The  optical  cross-section  is  then 


a 


(A15) 


In  the  more  general  case  where  (x,  y)  can  take  any  values,  the  above  double  integral 
can  be  efficiently  numerically  evaluated  by  first  approximating  the  integration  domain 
(that  is,  the  normalised  effective  aperture)  by  a  sequence  of  N  adjacent  trapezoids  [19]  as 
shown  in  Figure  A2.  Simple  analytic  expressions  for  the  above  diffraction  integral  over 
a  trapezoid  are  known  [20],  and  the  sloped  edges  of  the  trapezoids  better  approximate 
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the  edges  of  the  aperture,  £  =  ±F(rj)  as  the  number,  N,  of  trapezoids  is  increased.  The 
resulting  approximation  becomes 


N 

W(x,y )  «  !m(x,y) 

m= 1 


where  the  diffraction  integral  for  the  m-th  trapezoid  is 

—  {sin(£mz  +  rj my)  sine (%  +  awe]) 

+  sin(^mx  -  r)my)  sine  (h[y  -  amx ])}  x/0 

J-m\X,y)  = 

4 h  {£m  cos(r]my)  sine (hy) 

+otmrjm  sinc(2 r]my)  [cosh(y/2)  -  sinc(%)]}  x  =  0 


where 


sinc(a;) 

'Hm 

Ot-m 

h 


sin(o;/2) 

uj]2 


[1  w  =  0 

\  [F(mh )  +  F((m 

C m-\)h 

2  [F(mh)  —  F((m 
y/N. 


1  )h)] 

m] 


(A16) 


(A17) 


(A18) 

(A19) 

(A20) 

(A21) 

(A22) 


The  accuracy  of  the  above  approximation  improves  with  increasing  number,  N ,  of  trape¬ 
zoidal  divisions.  It  also  depends  on  the  shape  of  the  effective  aperture,  and  hence  on  the 
retroreflector  tilt.  Figure  A3  shows  plots  of  the  relative  error  of  the  numerical  approxima¬ 
tion,  calculated  over  a  range  of  values  of  N.  It  can  be  seen  that  the  relative  error  follows 
a  power-law  relation 

e  =  hN~k2  (A23) 

for  large  values  of  N,  where  the  values  of  the  parameters  k{  depend  on  the  tilt  angle.  In 
the  present  work,  N  was  chosen  to  be  1,500  in  order  to  give  values  correct  to  at  least  five 
significant  figures,  at  any  tilt  angle. 
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Figure  A3:  Relative  errors,  e,  in  the  numerical  approximation  calculation  of  the  optical 
cross  section  on  axis  ( 9r  =  0°),  for  which  an  explicit  expression  is  known  (equation  (A15)). 
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Appendix  B  Spoiled  cornercube  retroreflector  ray 

deflection 

Let  the  3-dimensional  coordinate  system  be  represented  by  the  three  orthonormal  vectors 
fiu  where  l  =  1,2,3.  Without  loss  of  generality,  the  three  planar  reflecting  faces  of  a 
cornercube  can  be  described  by  the  three  unit  vectors  that  are  each  normal  to  one  of  the 
faces,  and  by  fixing  the  apex  of  the  cornercube  at  the  origin.  The  three  unit  normals  can 
be  written 

Aj  =  (B1) 

i 

where  i  =  1,2,3  and  the  scalar  coefficients,  an,  are  the  Cartesian  coordinates.  In  the 
special  case  of  a  perfect  cornercube  whose  faces  are  all  mutually  orthogonal,  the  unit 
normal  vectors  are  equal  to  the  basis  vectors  so  that 

(Jil  =  Su  (B2) 

where  8n  is  the  Kronecker  delta  function.  More  generally,  for  an  imperfect  cornercube 

°u  =  A i  ■  A v  (B3) 

The  cosine  of  the  dihedral  angle,  0t],  between  adjacent  faces  whose  unit  normals  are 
and  Aj  wiH  be 


eij  —  A i  '  Aj 

=  'y  ,  O’il&jmp'l  ■  A  to 

(B4) 

=  VilGjl- 
l 


In  matrix  notation,  if  [e]  ^  =  etJ  and  \a]ij  =  aij  then 

e  =  aaT .  (B5) 


Note  that  the  diagonal  elements 


eu  =  l  =  Y^  (B6) 

i 

which  is  consistent  with  the  normalisation  of  the  /Gq.  A  perfect  cornercube  will  have 

eij  =  0  if  i  /  j. 

If  a  unit  vector  rh  =  (q  i  ,  0:9, 0:3)  represents  the  direction  of  a  ray  striking  surface  i  (that 
is,  the  face  with  unit  normal  /q)  then  the  direction  of  the  reflected  ray  will  be 

rhj  =  m  —  2(m  •  Ai)Ai  (B7) 

by  the  usual  law  of  reflection.  This  reflected  ray  will  in  turn  be  reflected  by  the  other  two 
faces.  The  order  in  which  the  ray  strikes  the  faces  will  depend  on  the  lateral  placement 
of  the  original  ray.  In  general  the  order  of  the  three  faces  can  be  denoted  i-j-k,  where 
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i  ^  j  /  k.  The  above  expression  for  the  single  reflection  can  be  repeatedly  applied  twice 
to  give  the  retro-reflected  ray  direction 

rhk  =  m-2^(rh- A{)Az 
i 

+  4  [eij(m  •  A'JAj  +  <Tfc(rh  •  Ai)A'k  +  ejk(*w  '  AjOAk] 

-8eij-eifc(m-Ai)Afc-  (B8) 

In  the  special  case  of  a  perfect  cornercube  the  terms  vanish  and  the  expression  sim¬ 
plifies  to  rhk  =  —  m,  as  expected.  For  an  imperfect  cornercube  the  terms  involving  a  dot 
product  can  be  written  in  terms  of  the  basis  vectors  by  using  equation  (Bl),  so  that 

(m  •  Ai)Aj  =  UlVilVjmVrm  (B9) 

l,m 


and  finally 


mk  =  rh-^2 

1,771 


2  ^  ^  Onl^nm  4  \€ij O’nO'jm  T  ^-ik^il^km  T  ^jk& jl&km\  T  ^^ij^jk^il^km  /  •  (BIO) 


The  deviation  of  the  retro-reflected  ray  from  the  perfect  case  is  then 

p  =  cos_1(— rh  •  rhk). 


(Bll) 


In  order  to  restrict  the  orientation  of  the  cornercube,  while  maintaining  the  three  degrees 
of  freedom  of  its  construction  that  are  its  dihedral  angles,  the  ot]  are  assumed  to  symmet¬ 
rical;  that  is,  (Tjj  =  aji .  If  the  normalisation  constraint  of  the  fi',  is  also  considered,  there 
remain  just  three  independent  components  of  the  a  matrix,  so  that  it  may  be  written 


^  ~  (?l+°b) 

V1  -  +  ac) 

y  crb  ac 


(?b 


\ 


\/l~  (ab+ac)j 


(B12) 


where  oa  =  rrl2/  ab  =  <713  and  ac  =  023.  It  follows  that  the  matrix  of  the  dihedral  angle 
cosines,  e,  is  also  symmetric  with  independent  elements 


ea  =  ei2  =  (Tay^l  -  (a%  +  of)  +  <Jay/l  -  (of  +  of)  +  a bac 

£b  =  £13  =  <*b  \/l  -  (of  +  of)  +  aacrc  +  abyJ  1  -  (of  +  of)  (B13) 

ec  =  e23  =  cracrb  +  <JcV  1  -  (of  +  of)  +  ac\J  1  -  (of  +  of) 

The  problem  at  hand  is  to  calculate  the  ray  deviation,  p,  for  any  given  values  of  the 

(ea,  eb,  ec).  To  do  this,  the  values  of  the  (aa,  ab,  <rc)  must  be  found  from  the  (ea,  eb,  ec)  using 
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the  nonlinear  equations  (B13).  This  can  be  done  numerically  using  any  of  a  number  of 
standard  nonlinear  root-finding  methods,  with  the  first-order  approximation 

(cra,  abl  crc)  ~  (ea,  eb ,  ec)/2  (B14) 

as  the  initial  estimate  of  the  solution.  The  value  of  the  ray  deviation  angle,  ip,  can  then  be 
calculated  using  equations  (Bll)  and  (BIO). 

A  further  consideration  in  calculating  the  retro-reflected  ray  direction  in  the  case  of  a 
solid  cornercube  is  the  refraction  at  the  front  surface,  both  on  entry  and  exit  of  the  ray. 
Let  the  front  surface  normal  be  fin  and  let  it  be  fixed  at 

An  =  (Al  +  A2  +  A3) /Vs.  (B15) 

The  incoming  ray  is  rh  =  a\jlx  +  02A2  +  03A3/  where  for  a  normal-incident  ray  a  1  = 
«2  =  03  =  —  1/ \/3.  The  ray  will  be  refracted  to  a  new  direction  m!  by  Snel's  law 

nm!  =  rh  —  Tpn  (B16) 

where  n  is  the  refractive  index  of  the  cornercube  material,  and 

T  =  y +  An  '  "I.  (B17) 

Similarly,  the  return  ray,  rh'k,  will  be  refracted  by  the  same  surface  to  give  the  ray  direction 

mk  =  nrh!k  +  T'  pn  (B18) 

where 

r'  =  \J  1  -  n2  +  n2(fin  ■  rh'k)2  -  n{fin  ■  rh'k).  (B19) 

Note  that  the  refracted  ray  directions,  m!  and  mk,  calculated  by  the  above  expressions 
are  not  normalised. 

It  is  now  possible  to  determine  numerically  the  value  of  the  dihedral  angle  cosine  ea, 
say,  as  a  function  of  eb  and  ec  that  gives  any  particular  values  of  the  retro-reflected  ray 
deviation  <p.  In  the  present  case,  the  ray  deviation  is  ip  =  3  arcsec  and  the  solid  cornercube 
material  N-BK7  has  a  refractive  index  n  =  1.519.  If  the  dihedral  angle  0„  differs  from  7r/2 
by  pa  then 

ea  =  cos  (7r/2  +  Pa)  (B20) 

with  similar  expressions  for  eb  and  ec.  By  considering  the  maximum  ray  deviation  pmax 
from  each  of  the  three  unique  permutations  of  the  reflecting  face  order  i-j-k,  the  two 
extreme  values  of  pa  as  a  function  of  the  other  dihedral  angle  errors  Bb  and  3C  can  be 
found  that  give  the  prescribed  maximum  ray  deviation.  The  results  of  these  calculations 
are  shown  plotted  in  Figure  Bl.  Here  a  normally  incident  ray  is  assumed.  It  can  be  seen 
that  in  the  present  case,  the  corresponding  tolerance  on  the  dihedral  angles  is  ±1.2  arcsec. 
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Figure  Bl:  Dihedral  angle  error  /3Q  as  a  function  of  the  other  dihedral  angle  errors,  /3b  and 
j3c,  that  gives  a  return  narrow-beam  deflection  equal  to  the  tolerance  (3  arcsec)  specified 
by  the  manufacturer  of  the  cornercube  retroreflector  considered  in  the  present  work.  The 
lower  and  upper  limits  are  shown  respectively  in  the  left-hand  and  right-hand  plots.  Val¬ 
ues  of  the  dihedral  angle  pair  (/3b,  /3C)  that  lie  outside  the  /3a  =  0  contour  lines  will  always 
exceed  the  ray-deviation  tolerance.  The  dihedral  angles  therefore  fall  within  approximately 
1.2  arcsec  of  the  (designed)  right  angles. 
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